	subroutine matvec(A,x,y,N)

	real*8 A(N,N),x(N),y(N)

	do i=1,N
		y(i) = 0.
	enddo

	do j=1,N
		do i=1,N
			y(i) = y(i) + A(i,j)*x(j)
		enddo
	enddo
	
	return
	end